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ABSTRACT 

The effects of gravitational darkening on the thermal structure of Be star 
disks of differing densities are systematically examined. Gravitational darkening 
is the decrease of the effective temperature near the equator and the correspond- 
ing increase near the poles of a star caused by rapid rotation. We also include the 
rotational distortion of the star using the Roche Model. Increasing the disk den- 
sity increases the optical depths in the equatorial plane, resulting in the formation 
'"q_j of an inner cool region near the equatorial plane of the disk. High rotation rates 

result in disks that have temperatures similar to those of a denser disk, namely 
cooler overall. However the effect of increasing rotation produces additional heat- 
ing in the upper disk due to the hotter stellar pole. Cool regions in the equatorial 
plane normally associated with high density are seen in low density models at 
high rotation rates. Gravitational darkening increases the amount of very cool 
and very hot material in the disk and decreases the amount of disk material 
at moderate temperatures. We also present models which study the effect of 
gravitational darkening on hydrostatically-converged disks, in which the temper- 
ature structure is consistent with vertical hydrostatic equilibrium. Because the 
equatorial regions become cooler, hydrostatically converged models that include 
gravity darkening have smaller vertical scale heights, and H/R is smaller by as 
much as 56 % near fcrit- Finally we explore differences in disk temperatures when 
alternate formulations of gravitational darkening, which lower the temperature 
difference between the pole and the equator, are used. 

Subject headings: Radiative transfer - stars: circumstellar matter - stars: emis- 
sion line. Be - stars: rotation 
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Introduction 



Jaschek et al.l (Il98l[ ) gives the current working definition of a Be star as "a nonsupergiant 



B star wfiose spectrum lias, or fiad at some time, one or more Balmer lines in emission." 
Tlie source of these emission hues is a geometrically thin, gaseous disk surrounding the 
star. These disks are often temporary, building up and dissipating over approximately a 
decade, although som e Be stars have shown persistent Balmer em ission over their whole 
observational history ( IWisniewski et al.l I2OIOI : ICarciofi et al.l |2012| ). When present, these 
disks produce the other characteristic fe atures of a Be star syst em. (1) A continuum excess 
occurs due to the cooler disk material ( jDougherty et al.lll99ll ). It begins as a weak, but 
noticeable, excess in the opti cal, resulting in con tributions to the (B-V) colour index of 
about a tenth of a magnitude (iHowells et al.l 1200 ll ). peaking at ~ lO/xm and continuing well 
past the infra-red (IWaters et al. The longest wavelength detection of a Be star is 6 cm 

( Dougherty et al. 199ll ). (2) Co ntinuum polarization is preseri t due to scattering in the non- 



spherical circumstellar material (IWaters fc Marlboroughlll992l ). (3) Emissi on occurs in other 



hydrogen line seri es, Lyman, Paschen, Brackett, Pfund, and Humphreys (iHony et al.l 12000 



Houck et al.ll2004l ). (4) Some Be stars show emissi on lines in elements other than hydrogen 
such as He I, Fe ll, and sometimes Si ll and Mg ll (jSlettebak et al.l Il992t iPorter fc Rivinius 
2003h . 



The formation mechanism(s) of these disks is still unclear, but it is likely facilitated by 
the rapid rotation of these stars which could allow a rel atively weak process to d rive mass-loss 
in a main sequence or just post-main sequence star (jPorter &: Riviniud l2003l ). Short term 
variability, on timescales between 0.5 to 2 c lays, is observed i n the photospheric line shapes 
and photometry of most classical Be stars (jPercy et al.lll994l ). and Be stars have also been 
known to exhi bit non-radial pulsations , indicating disturbances on and above the surfaces 
of these stars (jPorter fc Riviniud l2003l ). Disk building appears to a be stochastic p rocess 
occurring in a series of outbursts associated with increased emission (IStefl et al.ll201ll ). 



The B type stars that produce these disks are rapid rotators 



rotation rates are still a matter of debate ( ICranmerl l2005l : Ivan Belld |2012| ). but they are 



YudinI mm . Their 



fast enough that the effects of rotational distortion and gravitational darkening should be 
considered. As described in § 2 of our previous paper, iMcGill et al.l (1201 ll ). gravitational 
darkening is the rotationally induced reduction of the effective gravity towards the stellar 
equator causing a correspond i ng de creas e in the local temperature. The traditional formu- 
lation is f o und i n Ivon Zeipel I (119241 ) and ICoUind (119631 ). while newer results are discussed in 
van Belld (120121 ). Due to the small effect of rotation on the total luminosity, the pole of a 
rotating star is actually hotter than an equivalent, non-rotating star. Rapid rotation causes 
the star to become distorted with the radius becoming larger at the equator than the pole. 
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to a maximum of 1.5 times the polar radius. 



Interferometric observations of th ese stars have co nfirmed a variation in brightness across 
the surface of rapidly rotating stars ( van Bell3 2OI2I) . The polar regions of such stars are 



(van Bellel 


2Qm. 


(van Belle 


2012) 



noticeably brighter than the equator (Ivan Bellel |2012| ). However the temperature contrast 



(van Belle 


2012). 


0.25 ( 


von Zeipel 



Smaller values of indi cate a smaller 



temperature difference between the pole and the equator. I Che et al.l (120111 ) find their in 



terferometric observations of Regulus can be best explained by /3 = 0.188. Ob s ervat ionally 



determined values for /3 are typically between 0.25 and 0.18, and Ivan Bellel (|2012[ ) gives 
P = 0.21 as the typical value seen from interferometry. Interferometric data seems to sug- 
gest that this only af fects the tempe rature profile, while the shape of the star is consistent 



with a Roche Model (Ivan Bellel 120121 ). 



As the star is the energy so urce for the d i sk, ch anges introduced by rotation have an 
effect on the disk. As described in lMcGill et al.l (120111 ). the effects of gravitational darkening 
on models of classical Be stars is to reduce the temperatur e in the mid-plane of the disk 
while causing some temperature increases in the upper disk. iMcGill et al.l ( 120111 ) presented 
models for four different spectral types: BOV, B2V, B3V, and B5V at ten different rotation 
rates. However only a single disk density scale of po = 5.0 x 10^^^ g cm~^ (see Equation ([T|) ) 
was considered. In this paper, we examine the combined effects of variations in disk density 
and stellar rotation. We also include discussions on the effects of rotation on hydrostatically 
converged models and on the inclusion of different formulations of gravitational darkening. 
The organization of the paper is as follows: §|2] briefly outlines the models presented in this 
paper; § [3] provides our results; § 13.11 explores the effects of rotation when combined with 
changes in the density of the disk; § 13.2! looks at the effects of rotation on disk models that 
have been hydrostatically converged; § l3.3l examines the effects of using different formulations 
of gravity darkening on the disk temperatures; and conclusions are presented in § HI 



2. Calculations 



The modelling program, Bedisk, was used (ISigut fc Jonesll2007l ). Bedisk solves the 
statistical equilibrium equations for the atomic level populations and then enforces radiative 



ivon Zeipel I (|l924l ) defines (3 = 0.25 for radiative energy transport. Convective transport occurs at 
Tcs <^7000K which is less than 5 % of the total area our B5V model rotating at 0.99t;ciit so neglecting 
this should have no effect on our results. 
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Table 1: Adopted Stellar parameters 



Type 


Mass 


Polar Radius 


Luminosity 




l^crit 




"i^crit 




(Mo) 


(Re) 




(K) 


(rad 


) 


(km s^^) 


BOV 


17.5 


7.40 


3.98x10^ 


30000 


7.10x10- 


-5 


548 


B2V 


9.11 


5.33 


4.76 xlO^ 


20800 


8.38x10- 


-5 


466 


B3V 


7.60 


4.80 


2.58 xlO^ 


18800 


8.95x10- 


-5 


449 


B5V 


5.90 


3.90 


7.28 xlO^ 


15200 


1.08x10- 


-4 


439 



Notes: ujcnt 

^crit '"max 



8GM /27r^, where Vp is the polar radius of the star and M is its mass. 
Wcrit , where rmax 



1.5 Tp. Stellar parameters adopted from ICoxl (120001 1 



equilibrium at each poi nt of the computat ional grid representing the disk. The version 
of Bedisk described in iMcGill et al.l (l201l[ ) includes gravitational darkening and was run 
for the stellar parameters given in Table [T] and the rotation rates given in Tabl e [21 I n 



our ca lculations, we have assumed that the polar radius remains constant following I Collins 
(119661 ). The stellar temperatures are defined using the assumption that the total luminosity 
of the star rernains fi xed which is a reasonable approximation (for discussion of this point, see 
Lovekin et al.l ( 120061 )). Th ese stellar param eters were chosen to include a model from each 
of the five bins adopted by lCranmed (120051 ) to analyse the effects of spectral type on Be star 
rotational statistics. We have calculated sets of models for the spectral types BOV, B2V, B3V, 
and B5V. Unfortunately, B8V models are too cool for Bedisk when run for gravitational 
darkening near critical rotation without explicitly improving the treatment of the diffuse 
radiation field (i.e. disk self-heating) beyond the modified on-the-spot approximation used 



bv lSigut fc JonesI (l2007h 
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Table 2: Rotation rates. 



•^irac 


".'frac 


' eq/ ' p 




Vcq (km s ^) 










RflV 

-UU V 


dA V 


DO V 


UO V 


0000 


0000 


1.00 














0.2000 


0.2960 


1.01 


110 


93 


89.7 


87.7 


0.4000 


0.5680 


1.06 


219 


186 


179 


175 


0.6000 


0.7920 


1.14 


329 


280 


269 


263 


0.7000 


0.8789 


1.20 


384 


327 


314 


307 


0.8000 


0.9440 


1.27 


439 


373 


359 


351 


0.9000 


0.9855 


1.37 


493 


420 


404 


395 


0.9500 


0.9963 


1.43 


521 


443 


426 


417 


0.9900 


0.9999 


1.49 


544 


462 


445 


435 



Models presented in § 13.11 and § 13.31 assumed a fixed density structure for the disk in 
the form, 

p(i?,.)=pJ^Ve-(^/^(«))% (1) 



R 



where 



In Equation ([2]) for the disk scale height, fi is the mean- molecular weight of the gas (treated 
as a constant and set to half the mass of hydrogen) and Tiso is an assumed isothermal 
temperature used for the sole purpose of setting the density scale heigh t; Tjsn is set equa l 



to 60% of the effective temperature of the central star, listed in Table [T] (jSigut et al.ll2009l ). 
The form of the disk density given by Equation follows from the assumption of a radial 
power-law drop in the equatorial plane {z = 0) coupled with the assumption that the disk 
is in vertical hydrostatic equilibrium set by the z-component of the star's gravitational 
acceleration and the assumed temperature Tjso- This leads to a "fiaring disk" in which 
H oc 

Models were made with the density fall off par ameter, n, set to 3.0 . The disks were 



constructed with the unchanging grid described in iMcGill et al.l (120111 ) which keeps the 
disk empty for r < 1.5 rp. This allows a disk of fixed density structure to be used for all 
calculations as the star expands with rotation up to a maximum of 1.5 rp at critical rotation. 
The filled region begins at r = 1.6 rp and the density in the mid-plane of the disk at the 
beginning of this filled region is given by pbd = Po (1/1-6)^ = 0.24po using Equation ([T]). 
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The parameters po and pbd will be referred to as the density scale and the base density 
respectively. The models were constructed for four different values of the density scale, po'. 
5.0 X 10-1^ 2.5 X 10-1^ 1.0 X 10-1^ and 5.0 x lO-^^ g ^^-^ 



Models presented in § 13.21 have been hydrostatically converged which ensures that the 
vertical density structure is consistent with the temperatu re structure by enf orcing hydro- 
static equilibrium over an additional loop as described in ISigut et al.l ( l2009l ). Only eight 
hydrostatically converged models were made due to the long computation time required, 
10 times longer than an unconverged model. Using the BOV stellar parameters, two 
hydrostatically converged models were calculated for each density scale, {po'. 5.0 x 10^^^; 
2.5 X 10"^^; 1.0 X 10~^^; and 5.0 x 10"^^ g cm"^ ). The first model without rotation and 
the second, with a rotation rate of f = O.OQwcnt- 



Models presented in § l3.3l use a version of Bedisk that allows for variations in the formu- 
lation of gravitational darkening. Results for the B2V model were recalculated for the four 
fastest velocities in Table [2] at the four different density scales, po'- 5.0 x 10~^^; 2.5 x 10^^^; 
1.0 X 10^^^; and 5.0 x 10^^^ g cm^^. Two variants of gravitational darkening were used: 
(1) reducing the value of the exponent on the magnitude of gravity in von Zeipel's theorem, 
/3, from the standa rd value of . 25 to 0.18; (2) changin g the expression for the temperat ure 
structure based on Ivon Zeipel I (Il924j ) to that given by lEspinosa Lara fc RieutordI (l201ll ) . 
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3. Results 
3.1. Effects of Density Changes 




Fig. 1. — Cliange in four disk temperature diagnostics witli increasing stellar rotation and 
varying densities for the BOV model. Shown are the volume-weighted and the density- 
weighted average temperatures as defined in Equations ([3]) and (jl]) (upper left and upper right 
panels) and the maximum and minimum temperatures (lower left and lower right panels). 
Four different disk densities are shown with po = 5.0 x 10~^^; 1.0 x 10~^^; 2.5 x 10~^^; 
and 5.0 x 10~^^ g cm~^. Darker lines indicate higher po- The lower horizontal axis indicates 
the fractional rotational velocity at the equator and the upper horizontal axis indicates the 
corresponding fractional angular velocity. 
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Fig. 2. — Same as Figure [H but for the B2V model. 




Fig. 3. — Same as Figure [T] but for the B3V model. 
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Fig. 4. — Same as Figure [H but for the B5V model. 



Figures [H |2l [3l and H] show the changes produced by increasing rotation on the tem- 
perature of d i sks fo r a variety of densities using the same four temperature diagnostics as 



McGill et aLl (120111 ): the density- weighted average temperature, the volume- weighted aver- 
age temperature, the maximum temperature, and the minimum temperature. The density- 
weighted average temperature is defined as 



1 

and the volume-weighted average temperature is defined as 



Tp = -— I T{R,z)p{R,z)dV , (3) 

JWdisk 



Ty = -— T{R,z)dV. (4) 

l^disk J 

In order to avoid numerical effects, the median of the twenty hottest and twenty coolest disk 
locations are taken to represent the maximum and minimum disk temperatures. 

The most important trend seen in Figures [1] through H] is in the density- weighted average 
temperatures. There is a decrease in this temperature with both increasing rotation rate 
and increasing density. This trend is seen in all spectral types. In addition, for all spectral 
types, the curves for different po are well separated and decrease smoothly. For spectral 
types B2V, B3V, and B5V, the slope increases with increasing rotation but flattens out near 
critical rotation. For spectral types B2V and B3V, the four curves for each po are essentially 
parallel and for B5V are nearly parallel. Parallel curves are not seen in the BOV model, and 
the decrease in temperature with rotation is greatest for the largest density model. 
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The temperature minimums also decrease with both increasing rotation rate and in- 
creasing density, analogous to the density-weighted average temperatures. The size of the 
decrease in the minimum temperature is also larger for the least dense models. For spectral 
types B2V, B3V, and B5V, the curve for the the most dense model, po = 5.0 x 10~^^ g 
cm~^, is essentially fiat. This causes the temperature minimums to converge at fcrit to val- 
ues between 5000 and 6000 K for all spectral types. The curves for different po are distinct, 
but there is significant overlap. The temperature minimums in the non-rotating models are 
larger for the earlier spectral types, and the earlier models experience a larger decrease in 
temperature with both rotation and density. The drop in the value of the minimum tem- 
peratures with rotation and density is barely noticeable in the B5V models as all curves are 
nearly flat and cluster at psSOOO K. 

One of the most interesting effects produced by rotation on Be star disk temperatures is 
seen in the volume- weighted average temperatures of the B2V, B3V, and B5V models. These 
values initially decrease with moderate rotation but begin to increase approaching critical 
rotation. The strength of both the initial decease and the increase near critical rotation is 
larger for early spectral types, excluding BOV. For the two highest densities of the critically 
rotating B2V models, the volume-weighted average temperatures are the essentially the same 
as the non-rotating models. Generally the volume-weighted average temperatures are hotter 
for the least dense models, but the separations between the averages are not as large as 
for the density- weighted average temperatures or the temperature minimums, and there 
is a great deal of overlap. The effect of density is strongest in the B2V model and very 
small by B5V. For moderate rotation, lower density models generally have smaller volume- 
weighted temperature averages than denser models, but the lowest density B2V models 
(p = 5.0 X 10~^^), and the non-rotating and slowest rotating models of the second lowest 
density B2V models (p = 1.0 x 10^^^) do not follow this trend. For rapid rotation, lower 
density models generally have larger volume-weighted temperature averages than denser 
models, but this docs not occur for the B5V models. 

The volume-weighted temperature averages of the BOV model does not behave like the 
other spectral types. There is only a very small change in temperature with rotation, which 
is somewhat larger for larger densities. The least dense BOV model is of nearly constant 
temperature. The densest, non-rotating models are hotter than the least dense models, but 
by critical rotation, this trend has reversed. All of the values are tightly clustered between 
13900 to 14500 K. 

The maximum disk temperatures change very little with rotation or density. This is 
because these temperatures are found in parts of the disk that are not shielded by the dense, 
mid-plane and are directly illuminated by the pole of the star. The maximum temperatures 
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are determined essentially by spectral type. There is a small but noticeable decrease in 
temperature with rotation in all models, which flattens out near critical rotation. Like the 
volume-weighted average temperatures, the temperature maximums begin to increase for 
extreme rotation, significantly for the B3V models but very weakly for the other spectral 
types. 

One of the most complex behaviours seen in Figures [1] through H] is the how the re- 
lationship between the volume-weighted and the density-weighted average temperatures is 
affected by rotation and disk density. These effects can be separated into three categories 
based on density, the high and low density extremes and the case of intermediate density. 
(1) For the highest density BOV model, the two highest density B2V and B3V models, and 
the highest density B5V model, the volume-weighted average temperatures are always higher 
than the density-weighted average temperatures. This is due to the presence of a large, cool 
region in the equatorial plane at high density. (2) For the lowest density BOV models and 
the two lowest density B5V models, the volume-weighted average temperatures are always 
lower than the density-weighted average temperatures. This is due to the absence of a sig- 
nificant cool region in the equatorial plane at low density. (3) All models of intermediate 
density experience a transition: for low rotation rates there is either not a cool region in 
the equatorial plane or it is too small to cause the density-weighted average temperatures 
to be lower then the volume-weighted average temperatures. At higher rotation rates, grav- 
itational darkening causes the development of a larger cool region and the density-weighted 
average temperatures become less than the volume-weighted average temperatures. 

Figures O |6l [TJ and |8] show radial temperature profiles of these disks obtained by aver- 
aging the temperatures within each vertical column perpendicular to the mid-plane of the 
disk. The average temperature at each radial distance is density-weighted and is found using 

= rzl ' ■ (5) 

Jo p{R,z)dz 

As seen in the previous figures, lower density models have higher temperatures. The 
profiles follow two patterns: (1) for low density models, the temperature at the inner edge 
of the disk begins fairly hot, increases to a maximum and then the temperature decreases as 
the radius increases, sometimes reaching a constant value; (2) for sufficiently large densities, 
the temperature at the front of the disk still begins fairly hot, but the temperature drops to 
a minimum and then increases again at larger radii, sometimes reaching a maximum value 
and then dropping again. For all spectral types, the temperature minimum becomes cooler 
and the extent of the cool zone increases in size with rotation. The slope of the temperature 
increase at large radii becomes smaller with increasing rotation. If a temperature maximum 
occurs, its size is reduced with rotation and the slope of the temperature drop off at large 
radii becomes shallower. 
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Fig. 5. — Variation of the vertically- averaged, density- weighted average temperature with 
disk radius for various stellar rotation rates and disk densities for the BOV model. Each 
panel shows a different rotation rate: no rotation (upper left); v = O.SOfcrit (upper right); 
V = O.QOfcrit (lower left); and v = O.QQwcrit (lower right). Four different disk densities are 
shown for each rotation rate: with po = 5.0 x 10~^^; 1.0 x 10"^^; 2.5 x 10~^^; and 5.0 x 10~^^ 
g cm~^. The darker lines indicate higher po- 
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Fig. 6. — Same as Figure [5] for the B2V model. 





Fig. 8. — Same as Figure [5] for the B5V model. 



Figures [9], [TOl [TTl and [12] show vertical temperature profiles of these disks obtained by 
averaging the temperatures at different radii but at the same scale height, u = z/H{R). The 
average temperature is density-weighted and is found using 

J^""-- TiR,u)piR,u)AiR,u)dR 
j^-^^"- p{R,u)A{R,u)dR 
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where A{R,u) is the area function of the disk. A{R,u) is included in Equation [6] to account 
for the non-uniform spacing in r and z of the computational grid (i.e. Ar and Az both 
increase with R). The exact nature of this weighting matters less than the fact that all the 
models have been averaged in the same way. Comparing these plots illustrates the effects of 
rotation and density on the vertical structure of these disks. 

In these plots, the mid-plane of the disk is at u=0. All vertical temperature profiles show 
a temperature maximum occurring between one and two scale heights above the mid-plane. 
The location of this maximum moves higher above the mid-plane with both increasing density 
and increasing rotation. Rotation causes the width of the temperature peak to increase and 
become broader, while increasing density causes width of the peak to become narrower and 
sharper. Low density models without rotation have mid-plane regions either at the same 
temperature or hotter than the upper disk. Increasing the density causes the development 
of a cool region in the mid-plane. Rotation causes this cool region to form in the mid-plane 
at lower densities. At high rotation rates, pronounced cool regions are seen even in the least 
dense models. 

There is little variation in the temperature of the upper edge of the disk with either 
rotation or density in the B5V models and no significant change with density for the non- 
rotating BOV models. The upper edges of the disk increase in temperature with rapid 
rotation in the B2V and B3V models and have a minimum value at ffrac = 0.80. In the B2V, 
B3V and the rotating BOV models, the upper edges of the disk are hotter for lower densities. 
While there is a change in the temperature of the upper edges of the disk with density in 
the rotating BOV models, which increases in size with rotation, the median of this range is 
not significantly affected by rotation and remains at ~ 13000 K. 

Figure [12] shows the two dimensional temperature structure for a selection of eight 
BOV models. Two models are shown for each density, one without rotation and one for 
^^frac = 0.95fcrit- Clearly the temperature in the equatorial plane becomes cooler as the 
density increases; this has been noted many times. What is interesting is that this cool 
region appears at smaller densities when rotation is included. 
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Fig. 9. — Variation of the radially-averaged temperature profile for various stellar rotation 
rates and disk densities for the BOV model. Here the temperature is plotted versus the scale 
height, u, defined by u = z/H{R). Each panel shows a different rotation rate: no rotation 
(upper left); v = O.SOfcrit (upper right); v = O.QOfcrit (lower left); and v = 0.99fcrit (lower 
right). Four different disk densities are shown for each rotation rate with po = 5.0 x 10~^^; 
1.0 X 10~^^; 2.5 X 10~^^; and 5.0 x 10~^^ g cm~^. Darker lines indicate higher po- 
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Fig. 12. — Same as Figure 



|9]for the B5V model 
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Fig. 13. — Temperature profiles of tlie disk for a selection of BOV models. The left column 
shows models with no rotation. The right column shows models with a rotational speed of 
■^frac = 0.95. The top row shows models with po = 5.0 x 10^^^ g cm~^; the second row, 
Po = 1.0 X 10^^^ g cm^'^; the third row, po = 2.5 x 10^^^ g cm~^; and the bottom row, 
Po = 5.0 X 10^^^ g cm^^. The colour-map indicates the disk temperatures. The black line 
outlines the star. 
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Finally a useful way to illustrate the effect of density and rotation on the range of disk 
temperature is to construct histograms of T{R, z) for each model. Figures [H] and [15] each 
show a set of histograms for increasing rotation. Only two sets of histograms are shown 
for brevity, the B2V model with po = 5.0 x 10~^^ g cm~'^ (the lowest density considered) 
and another with po = 5.0 x 10~^^ g cm"'^ (the highest density considered). In Figure [HI 
the model without rotation is fairly warm and as rotation increases, a low temperature 
tail forms in the distribution. In addition, the fraction of disk temperatures in the highest 
temperature bins also increases at high rotation rates. In Figure the non-rotating model 
has a two strong peaks, one quite cool and a second middle peak. There is also a large and 
long high temperature tail and a weak high temperature peak because this model is dense 
enough to possess a cool region in the equatorial plane, with warmer regions above. As 
rotation increases, the middle temperature peak weakens and eventually disappears, while 
the low temperature peak becomes stronger and cooler. The weak, higher temperature 
peak becomes stronger. The distribution of temperatures in the fastest rotating model has 
two strong peaks at the temperature extremes. Overall, these histograms demonstrate that 
gravitational darkening increases the amount of very cool and very hot material in the disk 
and decreases the amount of disk material of intermediate temperatures. 
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Fig. 14. — Series of histograms of disk temperatures. All histograms plotted are for the B2V 
model with po = 5.0 x 10~^^ g cm~^. Each panel shows a model with a different velocity as 
indicated in the left of each panel. The upper most panel is the non-rotating model and the 
rotation rates increase downward. The lower most panel is the fastest rotating model with 
■^frac = 0.99fcrit- The filled circles indicate the density-weighted average temperatures and 
the filled squares indicate the volume-weighted average temperatures for each model. 
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Fig. 15. — Same as Figure [T4l but with po = 5.0 x 10~^^ g cm~^. 



It is often useful in the study of circumstellar material to compare the temperatures 
found in the disk to the effective temperature of the star. However, when the star is a rapid 
rotator, there is a range of surface temperatures, and some representative value must be 
found. The most basic definition of effective temperature is from the Stephan-Boltzmann 
law, namely T^g = (L/aA)^^^ where A is the surface area of the star. Extending this to a 
rotating star gives, Tefr(ffrac) = {L/aA{vbac))^^'^- If the luminosity is considered unaffected 
by rotation we have. 



^efr(^^frac) = Tcs{V{re,c = 0) 



^(^frac = 0) 
^(ffrac) 



1/4 



(7) 



Figure [16] shows the density-weighted average temperatures divided by the stellar effective 
temperatures defined by Equation ([7]) versus rotation rate for all models. This plot summa- 
rizes the predicted global temperatures of the disks around rapidly rotating stars of a wide 
range of disk densities. The temperature ratio of the disk to the star is between 0.40 - 0.65 
for all models. Denser disks are proportionally cooler. The disks around cooler stars all have 
temperatures which are larger fractions the stellar effective temperatures than those around 
hot stars. The ratio between the density-weighted average temperatures in the disk and the 
effective stellar temperatures drops with moderate rotation. Approach to critical rotation 
causes an increase in the ratios with rotation in both cooler stars and less dense disks. The 
slope of the ratio flattens out for rapid rotation in hot and dense disks. In all cases, the 
trends are dominated by the variations in disk density, not rotation. 
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Fig. 16. — Range of density- weighted average temperatures for the all model disks divided by 
a measure of stellar effective temperatures including the effects of rotation (see Equation ([7])) 
versus rotation. The median temperature for the four different densities is plotted. The 
upper error bar corresponds to the highest value (and the lowest density model) and the 
lower error bar corresponds to the lowest value (and the highest density model). Each panel 
is a spectral type as indicated. The lower horizontal axis indicates the fractional rotational 
velocity at the equator and the upper horizontal axis indicates the corresponding fractional 
angular velocity. 



Table [3] gives the density- weighted temperatures for all density scales and spectral types 
without rotation and at maximum rotation {v = 0.99fcrit)- This table shows the full range 
of temperature changes that occur with rotation for different densities. Increases in either 
density or rotation cause decreases in the density-weighted temperatures. For the BOV 
models, these two effects strengthen each other; denser disks are more strongly affected by 
increasing rotation and disks surrounding critically rotating stars are more strongly affected 
by density changes than disks surrounding non-rotating stars. The density-weighted average 
temperatures change by —6 % as rotation varies from zero to 0.99ucrit for the least dense 
models, and change by —21 % for the densest model over the same range in velocity. When 
we look at the results with rotation rates fixed, we find that the density-weighted average 
temperatures change by —14 % from = 5.0 x 10^^^ g cm~^ to po = 5.0 x 10^^^ g cm~^ 
without rotation, and change by —28 % for the fastest rotating models. In the B2V and 
B3V models, the effect of rotation is nearly the same for all po, and the percent differences 
range between —13 and —15 % for the B2V models and between —15 and —16 % for the 
B3V models. The effect of density changes is also nearly same on models without rotation 
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and those with a rotation rate of 0.99T;crit, with the percent differences ranging between —23 
and —26 % for the B2V models and remaining unchanged at —24 % for the B3V models. 
By B5V the these two effects weaken each other. The least dense models are more strongly 
affected by rotation and the non-rotating models are most affected by density changes. The 
density- weighted average temperatures change by —20 % from zero to 0.99fcrit for the least 
dense models, and change by —12 % for the densest model. Holding rotation rates fixed, we 
find that the density- weighted average temperatures change by —24 % from po = 5.0 x 10~^^ 
g cm~^ to Po — 5.0 X 10~^^ g cm~^ without rotation, and change by —16 % for the fastest 
rotating models. 

Table 3: Density-weighted awra.go tomporaturos without rotation and with Vf^ac — 0.99 



Type 


Po 


2;(o.oo) 




% A(Vfra^) 




g cm~^ 


K 


K 




BOV 


5.0 X 10-12 


15200. 


14250. 


- 6. 




1.0 X 10-11 


15180. 


13920. 


- 9. 




2.5 X 10-11 


14520. 


13010. 


-11. 




5.0 X 10-11 


13280. 


10780. 


-21. 




% A(p) 


-14. 


-28. 




B2V 


5.0 X 10-12 


11620. 


10220. 


-13. 




1.0 X 10-11 


11090. 


9690. 


-13. 




2.5 X 10-11 


10200. 


8800. 


-15. 




5.0 X 10-11 


9190. 


7870. 


- 15. 




% A(p) 


-23. 


-26. 




B3V 


5.0 X 10-12 


11000. 


9520. 


-15. 




1.0 X 10-11 


10590. 


9110. 


-15. 




2.5 X 10-11 


9670. 


8230. 


-16. 




5.0 X 10^11 


8690. 


7460. 


-15. 




% A(p) 


-24. 


-24. 




B5V 


5.0 X 10-12 


9710. 


7970. 


-20. 




1.0 X 10-11 


9290. 


7760. 


-18. 




2.5 X 10-11 


8450. 


7230. 


-16. 




5.0 X 10-11 


7660. 


6760. 


-12. 




% A(p) 


-24. 


-16. 





Taking all the results into consideration, we see that while the bulk of the disk becomes 
cooler, there is evidence of heating in some parts of the disk (either from the volume- weighted 
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average temperatures or the maximum temperatures or both). Heating of the upper edge 
of the disk is seen in the vertical profiles of the B2V and the B3V models, as shown i s 
Figures [TUl and [TTl This was also seen in Figures 8, 9 and 13 - 15 in iMcGill et al.l ( 1201 ll ). 
When density changes are combined with rotation, we see that the cooling associated with 
increasing density occurs at lower po at higher rotation rates. This suggests that if rotation 
is not taken into account, the densities of these disks could be over estimated. 



3.2. Self-Consistent Vertical Hydrostatic Equilibrium 




Fig. 17. — Ratio of disk scale height to disk radius versus disk radius for the fixed density 
structure shown in Equation ([1]), and two hydrostatically converged models, one without 
rotation (ffrac = 0) and one with rotation including gravitational darkening (ffrac = 0.99) 
for the BOV model. Four different disk densities are shown: po = 5.0 x 10~^^ (upper left); 
1.0 X 10"^^ (upper right); 2.5 x 10^^^ (lower left); and 5.0 x 10"^^ g cm^^ (lower right). 



In this section we describe the effe cts of rotation on disk models that have been hy- 
drostatically converged, as described in ISigut et al.l ( 120091 ). The vertical structure of clas- 
sical Be star disks are believed to be in hydrostatic equilibrium. Equation ([1]) governs the 
density structure of a hydrostatically supported isothermal disk, but real Be star disks are 
not isothermal. Temperatures throughout the disks typically vary by factors of 2 to 3, as 
seen in Figures [H through H] and Figures [TH and [IB This means that Equations and ([2]) 
are inconsistent with the detailed T{R, Z) distributions for the models computed from ra- 
diative equilibrium. Vertical pressure support in the disk requires that the vertical density 
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profile follow —dP/dz = p{z,r)gz, in which Qz is the vertical component of the star's gravi- 
tational acceleration. The local pressure is strongly affected by the local temperature of the 
disk via the equation of state, P = pkT/fi. This means that in order to create a density 
profile consistent with the temperature solution, an additional loop is required within the 



Bedis k code to enforce hydrostatic equilibrium in each column, as described by ISigut et al. 



(120091). Because of the extra loop, a converged model takes approximately ten times longer 



to calculate than an unconverged model. 

Figure [T71 shows the ratio of the disk scale height to the disk radius, H{R)/R, versus disk 
radius, R, in the inner disk (r < 5rp) for BOV models of both the fixed density structure 
described by Equation ([1]) and those that have been hydrostatically converged. H{R) is 
the height where the density drops by a factor of 1/e. For an isothermal disk, H{R)/R is 
equal to the ratio of the sound speed, Csound = \/ Pj P-, to the Keplerian orbital speed, forbit = 
\l GMjR. Because of the nature of the fixed isothermal density structure, it always produces 
the same if profile that increases as R}!"^ regardless of po or n\ as it is a function of only 

Tiso- In the current set of models, H{R)/R is always larger for the fixed structure compared 
to the hydrostatically converged structure because we have adopted the usual result that 
^iso = O.GTeff. However this temperature is often higher than the temperatures actually 
found in the mid-plane of den se disks. Th e refore , Equation ([1]) over-estimates the pressure 



support of the mid-plane, (see ISigut et al.l (120091 ) for more details). The difference between 



Tiso and the actual disk temperature becomes smaller as the disk radius increases and this 
is why the difference in H{R)/R between the hydrostatically converged models and those 
with the fixed density structure decreases with increasing radius. The difference in H{R)/R 
between the fixed structure and the hydrostatically converged models also increases as po 
increases because the temperatures in the mid-plane decrease strongly as po increases. Note 
that there is no single choice for Tiso that would be able to match the correct vertical density 
structures of the hydrostatically converged models shown in Figure [T71 

We now demonstrate that rotation, in addition to density, affects the scale heights 
computed for these disks. Figure [17] demonstrates that H{R)/R for v = 0.99 ffrac is always 
lower than that for v = for the same radius and density. This is because rapidly rotating 
models have cooler temperatures in the mid-plane, resulting in smaller scale heights: rapid 
rotation causes models with thinner disks. The model with po = 5.0 x 10~^^ g cm~^ and 
f = 0.99 f crit has a very thin disk with a maximum value in the inner disk of H{R) /R = 0.033 
at 2.6 stellar radii and a minimum of 0.027 at 3.75 stellar radii. 

Table m shows the density- weighted average temperatures without rotation and with ro- 
tation at 0.99vfrac for the fixed density structure and the hydrostatically converged structure 
for the four different densities shown in Figure [T71 Increasing density or rotation causes a de- 
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crease in the density-weighted average temperature with or without hydrostatic equihbrium. 
Requiring self-consistent hydrostatic equihbrium reduces the effects of changing rotation and 
density and also increases the temperatures found in the dense BOV models. 



Table 4: Density-weighted average temperatures (K) with and without self-consistent hydro- 
static equilibrium 



BOV 




Fixed 






Hydro 




Po (g cm"^) 


^^frac = 0. 


^^frac = 0.99 


% A(t;frac) 


^^frac = 0. 


Vlra.c = 0.99 


% A(t;frac) 


5.0 X 10-12 


15200 


14250 


-6 


15170 


14300 


-6 


1.0 X 10-11 


15190 


13920 


-9 


15120 


14070 


-7 


2.5 X 10-11 


14520 


13010 


-11 


14590 


13500 


-8 


5.0 X 10-11 


13280 


10780 


-21 


14020 


12570 


-11 


% A(p) 


-14 


-28 




-8 


-13 





3.3. Different forms of gravitational darkening 



As described in Ivan Belld ( 12012| ). gravitational darkening has been interferometrically 
confirmed, but the difference in brightness over the stellar disk is not as strong as tha t 
predicted by the standard formulation of gravitational darkening presented by I Collins! ( 119631 ). 
This formulation states that 

F = C^g'^, (8) 

where F is the local radiative flux and C^j is constant across a star and determined by the 
luminosity such that /. 

C^ = LJ g^^dA, (9) 



where the integration is over the surface of the star. The luminosity, L^, caru b e tr e at e d 
as a c onstant or as a function of rotation. The canonical value of /3, deflned in Ivon Zeipel 



(Il924l ) , is 1/4, but ob s ervations sugg est a range of values between 0.25 and 0.19 for B type 
stars (jvan Belle 2012 ). Claret ( 2012 ) presents a calculation for the variation of /3 with local 
temperature and optical depth in a rotating star of 4 Mq and describes departures from 
von Zeipel in the upper layers where f3 can be as low as ~ 0.18 when energy transport is 
radiative. Figure [TS] shows the variation in temperature with rotation for P = 0.18. and 0.25 
at four different stellar co-latitudes. The constant C^j was recomputed in each case. Using 
P = 0.18 weakens both the temperature increase in the polar region and the temperature 
decrease in the equatorial regions, reducing temperature difference between the equator and 
the pole. 
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Fig. 18. — Surface temperature versus rotation rate for four different co-latitudes, 0", 30°, 
60", and 90°, for the B2V model. The solid lines indicate (3 = 0.25 while dots indicate 
/5 = 0.18. The lower horizontal axis indicates the fractional rotational velocity at the equator 
and the upper horizontal axis indicates the corresponding fractional angular velocity. 
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Fig. 19. — Change in disk temperature diagnostics with increasing rotation for the B2V 
model with po = 5.0 x 10~^^ g cm~^. The lines marked with upward pointing triangles indi- 
cate the maximum temperature and the lines marked with downward pointing triangles indi- 
cate the minimum temperature. The lines marked with circles indicate the volume-weighted 
average temperature and the lines marked with squares indicate the density- weighted average 
temperatures. The black and grey lines show models with [3 = 0.25 and ,6 = 0.18 respec- 
tively. The lower horizontal axis indicates the fractional rotational velocity at the equator 
and the upper horizontal axis indicates the corresponding fractional angular velocity. 
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We choose /3 = 0.18 as a lower limit and re-ran Bedisk for the four highest rotation 
rates assuming the B2V star surrounded by a disk with po = 5.0 x 10~^^ g cm~^. The results 
are shown in Figure Uni Lowering /3 reduces the strength of gravitational darkening and the 
corresponding effect of the stellar flux on the temperatures in the disk. The density-weighted 
average temperature decreases with rotation and this decrease is less with f3 = 0.18. The 
volume-weighted average temperatures in the disk have a more complex behaviour, initially 
decreasing then increasing with rapid rotation, as /3 is reduced to 0.18. The initial drop in 
temperature becomes smaller and the increase in temperature at higher rotation rates is also 
smaller. The value of /3 does not effect the maximum or minimum temperatures found in 
the disk. 

The discrepancy in /3 between theory and observation has motivated a closer look at 
gravit ational darkening and possible re finements to the theory. lEspinosa Lara fc Rieutord 
( 20111 ) presents such a refinement. Both lEspinosa Lara fc RieutordI ( 1201 ll ) and the canonical 
method begin with the expression for radiative energy transport. The equation of hydro- 
static equilibrium is used to relate the radiative fiux to the local gravity, resulting in ex- 
pres sions of the form, Fr^A oc q, whic h becomes Teg oc \g\^^'^ using the Stephan-Boltzmann 
law ( von Zeipel 1924 : Clayton 19831). The difference in these approaches is how the other 
terms in the equation are handled. Ivon Zeipel I (Il924[ ) treats the equipotentials as isobaric 
surfaces which implies that all other terms in the equat ion are c onstan t over the surface 
of a rotating star (for this type of argument see §6.8 of IClaytonl (jl983[ )). This leads to a 
contradiction: the gas temperature is taken as a constant across the surface of a rotating 
star, while the effective temperature decreases with increasing stellar co-latitude. Although 
this is not truly a paradox (as the effective temperature and the gas temperature are not the 
same quantities), it would be reasonable for a higher radiative flux at the poles to produce 
higher local temperatures. 



Alternatively, lEspinosa Lara fc RieutordI (l201l[ ) make no assumptions about the other 
terms necessarily being constant, letting F^ad = f{r,(^)g- The unknown function, f{r,6), is 
found by requiring that V--Frad = 0. Using the Roche model and solving the resulting partial 
differential equation gives 

. L,. 



1/4 



eff 



aGM' 

where 9w is deflned by the requirement that 



tan 9u 
tan 6* 



-9 



1/4 



COS Bon + In tan — = —x w cos 6 + cos 6 + \n tan - 
2 3 2 



(10) 



(11) 



In Equation [m 9 is the spherical coordinate, x is a scaled radius r/req(ix'star)! and w is a dif- 
ferent fractional angular velocity given by u^tax \ rl^{bJs\,ai)/GM which is notutnc- Equations 
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of this form which describe a variable f{6) are sometimes written in the equivalent form of 
T = foQ^"'^^^^-' , where fo is now constant with all variation accounted for in S(9). Gravitv 
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Fig. 20. — Ratio of the temperatures predicted bv lEspinosa Lara fc RieutordI ( 2011 '). Telr 



to those of standard gravitational darkening, Tyz, ( Ivon Zeipel I Il924j : ICoUinsI Il963l ) versus 
rotation rate for four different co-latitudes, 0°, 30°, 60", and 90". The lower horizontal 
axis indicates the fractional rotational velocity at the equator and the upper horizontal axis 
indicates the corresponding fractional angular velocity. For readability we only show rotation 
rates faster than 0.75wfrac, since this ratio is ~ 1 for low rotation rates. 



Equation (fTO)) pre dicts lower polar temperatures and higher equatorial temperatures 

than von Zeipel I ( 1924 ). This is similar to, but not the same as, reducin g (3. At slow rotation 

rates t here i s very little d i fferen ce between the temperatures predicted by lEspinosa Lara fc Rieutord 



(I2OIII ) and Ivon Zeipel Jl924h . but the differences increase with rotation and the ratio 



between them actually diverges at the equator for critical rotation because while both 
functions find a temperature of zero at the equator for critical rotation, the prediction of 
Espinosa Lara fc RieutordI (120111 ) goes to zero with rotation more slowly (see Figure [20]) . Fig- 
ur e [2T] shows the effe ct of using th e Equa tion (fTOj) compared to traditional gravity darkening 
of Ivon Zeipel I ( 1l924j ) and ICoUinsI ( 1l963l ) . It is similar to Figure [191 predicting a weakening 
of the effects of gravitational darkening. However, Equation (fTOj) reduces the heating of the 
pole more than the cooling of the equator, in comparison to simply reducing /5. This is why 
using Equation (ITU!) does not effect the density-weighted average temperatures as much as 
setting P = 0.18. Using Equation (|TOl) has a similar effect on the volume-weighted average 
temperatures as setting (3 = 0.18, reducing the size of the initial drop in temperatures as well 
as reducing the heating at high rotation rates but the effect of setting /3 to 0.18 is stronger. 
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Fig. 21. — Same as Figure [T9] except the black lines show standard gravitationa l darke ning 
with (3 = 0.25 while the grey lines show results using lEspinosa Lara fc RieutordI ( l201ll ). 



4. Conclusions 

Both density and rotation significantly affect disk temperatures. This is clearly seen 
in the density-weighted average temperatures in Figures [1] through H] and in Table [31 Both 
cause disks to become cooler. Density is a stronger controller of the temperatures of these 
disks than rotation. This is as expected. Increasing the optical depth through the disk 
decreases the amount of photoionization radiation able to penetrate the disk by a large 
factor, especially in the mid-plane where the bulk of the material is located. Alternatively, 
rotation does not change the luminosity of the star, but only redirects it away from the 
mid-plane. The effects of rotation on the thermal structure of disks are not noticeable below 
0.20 f frac- For moderate rotation, 0.20 f crit to 0.60 fcrit, there are small but noticeable changes 
in the disk temperatures. From 0.60fcrit to O.SOwcrit the effects of rotation become stronger 
as the equator cools and from O.SOfcrit to 0.99fcrit the stellar pole becomes hot enough to 
influence the disk while the equatorial cooling continues. 

Increasing rotation from zero to O.SOffrac can have the same effect on the density- 
weighted average temperatures as increasing the density by a factors of 1.5 to 5, depending 
on the model. The effects become even larger closer to critical rotation. Increasing rotation 
from zero to 0.99ffrac can have the same effect on the density-weighted average temperatures 
as increasing the density by 2.5 to 7.5 times, depending on the model. Therefore, while not 
as strong a temperature controller as density, rotation can be very signiflcant for moderate 
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to strong rotation and should not be neglected. Be cause classical B e stars rotate faster than 



0.40ffrac; ^ud many rotate faster than O.SOffrac (jCranmerl |2005| ). gravitational darkening 



should be included in models of classical Be stars disks . 

With increasing rotation and density, Be star disks become less isothermal. This is can 
be clearly seen in Figures [H] and [151 These histograms show an increase in the amount of 
very cool gas as rotation increases, and at high rotation rates, there is an increase in the 
amount of hot gas also. Increasing the density of these disks increases the amount of cool 
gas, but does not significantly increase the temperatures in the upper disk. Therefore, only 
strong rotation causes an noticeable increase in the temperatures of the upper disk due to 
the hotter stellar pole. 



Classical Be stars are known to be rapid rotators with geometrically thin disks ( iPorter &: Rivinius 



20031 ). The inclusion of gravitational darkening in hydrostatically converged models shows 
that these disks are predicted to be very thin around rapidly rotating stars because the large 
equatorial cool region reduces pressure support in the vertical direction. The scale heights 
predicted for the densest model considered are very small, with H/R reaching as low as 0.027 
at 3.75 stellar radii. 

Unsurprising, using weaker forms of gravitational darkening weakens the effects of gravi- 
tational darkening on the disk. Simply reducing /3 decreases the effects of rotation on both the 
volume-weighted average temperatures and the density-weighted average temperatures in the 



disk, b ut also introduces a new free parameter. The formulation of lEspinosa Lara &: Rieutord 



(120111 ) is a excellent alternative as it offers a physical explanation, avoids adding an artificial 
parameter, and is an algebraic solution so that it can be utilized for any star and any rotation 
ra te lower than critical. Rep l acing standard gravitational darkening with the formulation 



of lEspinosa Lara fc RieutordI (120111 ) does not change the density- weighted average temper- 
atures in the disk very much, but it does reduce the heating of the upper disk due to the 
stellar pole. There are other physical phenomena which could further effect the te mperature 



struct ure of a rotating star, such as horizontal transport of energy as described in iHadrava 



(I1992I ) which would re-distribute energy from the bright poles to the dimmer equator, but 



this effect is expected to be small. 

In the future, we plan to produce observables for our rotating Be star disk models. It 
will be interesting to see how rotation and density changes effect the colours and Ha emission 
from these systems, which are often used as the primary diagnostics of our circumstellar disk. 
For disks with large optical depths, temperature changes can strongly influence disk emission 
so effects due to rotation are likely to be important. 
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